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Abstract 

Vibrational properties of ultralong Rydberg Cs2 molecules are investigated on corresponding potential 
curves obtained by perturbation theory. The Rydberg Cs2 molecules are associated by a Rydberg Cs(nS/nP) 
atom (n = 30 — 70) and a ground state Cs(6s) atom. The starting point for the perturbation treatment of cor- 
responding Rydberg molecular potential curves is to generate accurate atomic Rydberg states from realistic 
ab initio effective core potential. The calculated results have similar characteristics with available exper- 
imental and theoretical investigations on Rydberg Rb2 molecules. And this is the first time that Rydberg 
molecules are studied at the ab initio level. 
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I. INTRODUCTION 



Owing to the rapid development on ultracold atoms and molecules [l 



-bi, 



a novel type of 



ients| J, 0]. 

molecule is typically formed by the unusually long-range interaction between a Rydberg atom and 
a ground state (or another Rydberg) atom. Consequently Rydberg molecules are classified into 
two types according to their formation. In the first case, a macrodimer Rydberg molecule Jf]-^! 
associated by two Rydberg atoms (A*+B*) has a size larger than 1 /jm and in the second case a 
butterfly or trilobite like Rydberg molecule 111 Oh 1311 associated by a Rydberg atom and a ground 
state one (A*+B) has a typical size of about 10 3 ao (Bohr radius). Since the first calculation of 
Rydberg molecular potential curves in 2000 by Greene and coworkers ll 1011. Rydberg molecules 



141—1 1 VM due to their ex- 



have been extensively investigated both theoretically and experimentally [ 
traordinary properties such as huge size, large dipole moment, weak bounding energy and long 
lifetime. In particular, in some fast developing area of quantum physics, Rydberg molecules are 
more suitable than Rydberg atoms for quantum manipulation with electric fieldfl la. 1 19 h 

The A*+B* and A*+B types of Rydberg molecules have been postulated in Refs. |7] and Jioll 



respectively. Soon experimental evidence Il20rl22ll has been found for the A*+B* type of Rydberg 



23] 



molecules. Later on experimental verification of the A*+B type of Rydberg molecules fl 14 
has been reported. In 2009 both types of Rydberg molecules have been reported to be directly 
observed experimentally fllsj]. In Ref. fl Bendkowsky and coworkers have reported the measured 
spectra of Rydberg molecular states formed by Rb(5s) and Rb(Hs) with high principle quantum 
number n = 34, • • • ,40. While in Ref. Q5[] electric field induced Cs Rydberg atom macrodimer s 
bound at internuclear separations of 3 — 9 /urn are investigated, focusing on (63d) + (65 d), (64d) + 
(66d), (65 d) + (61 d) and (66d) + (68 d) pairs. Soon after that Rydberg trimers lll5ll have also 
been created by a single-step photoassociation. For a better understanding of the mechanism 
and further manipulations of Rydberg molecules, coherent transfer of initially free pairs of Rb(5^) 
atoms to Rydberg molecules (Rb*+Rb) have been demonstrated experimentally Q24J] . Very recently 
permanent electric dipole moment of a homonuclear Rydberg molecule has been observed ||18|l to 
has the order of 1 Debye even for highly symmetric case of Rb(n5)+Rb(55) which further pave the 
way for possible experimental manipulation of this type of molecules. 

From theoretical point of view, there have been some calculations of the interaction potential 



curves and vibrational properties of Rydberg molecules [|7 



Jul In, 



1911 most of which focusing on 
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A*+B* type. The available theoretical calculations have played significant roles in the experimen 



tal realization and interpretation of molecular Rydberg states. In Ref. Ill7h the photoassociation 



rates of different bound states have been investigated for Rb Rydberg microdimers which implies 
potential application to quantum information processing. Very recently the electric field control of 
triatomic Rydberg molecules have been reported Jl9ll which shows their sensitivity to external elec- 
tric field and non-adiabatic transitions between neighboring Rydberg molecular potential curves. 
For A*+B type of Rydberg molecules, in general the available theoretical investigations Jiol 11] 



have adopted s wave scattering or perturbation treatment based on Rydberg states from certain 
empirical electron-core potential for the Rydberg atom. Therefore accurate description of the Ry- 
dberg atom is essential for the accuracy of the Rydberg molecular potential curve no matter by 
scattering or perturbation treatments, in particular when pursuing quantitative agreement with ex- 
periments. To meet this challenge, the present work will adopt combined ab initio calculations 
and a long range analytical form for better description of the electron-core potential of a Rydberg 
atom to obtain accurate potential energy curves and vibrational properties of a Rydberg molecule. 
According to our best knowledge this is the first time that Rydberg molecules are studied at the ab 
initio level. 

In this work ultralong Rydberg CS2 molecular potential curves dissociating into 
Cs(nS / nP)+Cs(6s) and corresponding vibrational levels will be investigated (n = 30 — 70). De- 
tails of the methods including combined ab initio and analytical description of Cs Rydberg atoms, 
perturbation treatment of Cs2 molecular potential curves and diagonalization of Cs2 vibrational 



Hamiltonian in Fourier grid discrete variable representation [|25|l are given in Sections Til Al III Bl and 
III CI respectively. Then the Rydberg Cs2 molecular potential curves, vibrational levels, equilib- 
rium internuclear distances and ground state binding energies are discussed in Section [TTTJ Finally 
the conclusions of the present work are summarized in Section ITVl 

II. METHODS 

A. Hamiltonian and Wavefunction of Cs Atom 

The probability density of Rydberg electron is crucial for the potentials of the long-range Ry- 
dberg molecules. Therefore the Rydberg states of a cesium atom need to be accurately solved 
as a prerequisite for studying Rydberg Cs2 molecules. Based on central-field approximation, the 
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cesium atom can be described by a single electron model with the core-electron interaction using 
an effective core potential. The Hamiltonian of a cesium atom reads 

h 1 

#atom = -^V 2 + V atom (r), (1) 

2m 

where V a tom(/) is the effective core potential. For large r the effective core potential has standard 
long range asymptotic behavior as 

e e a c 

4n£ r (4tc£o) 2 2r 4 

with a c being the Cs + core polarizability. While for small r the potential becomes considerably 
more complicated consequently the short range potentials are directly obtained by ab initio calcu- 
lations using Gaussian09 suit of pro grams . Il26ll The CCSD(T) method, coupled cluster including 
single and double excitations and triple excitations by perturbation, is adopted to get accurate 
energies. For modern quantum chemistry calculations the CCSD(T) method has been widely 



used for accurate energy calculations of small molecules[|27|,|28|] due to its balance between speed 
and accuracy. In this work all ab initio energy calculations are performed with Gaussian09 at 
CCSD(T)/QZVP level of theory and full correlation between all the electrons explicitly involved 
in the QZVP basis set has been calculated. 

Having the Cs atom Hamiltonian at hand, the single electron Schrodinger equation can be 
solved to get the valence electron wavef unctions. The total wavefunction ^atomtr) can be 
separated into the radial and angular parts as l I , atom(r) = ^atomi/, 6, cp) = JHn/(r)I^ OT (9,<p) with 
F, m (8,<p) being the well known spherical harmonics. While the radial wavefunction %K n i{r) can 
be solved by the numerical methods such as Numerov method 12911 • Taking X(r) = r£H /I /(r), the 
radial Schrodinger equation can be written as standard Numerov form = —f(r)X, where 
f(r) = [Eni — Vatom(f)] — ^jr^ w ^h the corresponding Cs atom eigenenergy E n \. From the 
quantum defect theory, the fine structure eigenenergy of Cs atom E n ij can be obtained with the 
quantum defects reported in Ref. 113011 . With E„ij it is straightforward to derive the eigenenergy 
E n [. Discretizing r by r,- = (N — j + 1) • h, where h is the step size and N is the number of total 
grids, the Numerov integration is written as 

x i+1 = ± 6 -J- A — — I + (h 6 ). (3) 

1 + 13/1+1 

The final integration is being done starting from the long range limit where the eigenfunction has 
asymptotic formXi ? = exp (— ^/ —2mE n ir 1^2/ h) . 



B. Perturbation Treatment of Rydberg CS2 Molecular Potential Energy Curves 



A Rydberg Cs2 molecule, formed by a ground state Cs(6s) atom and a Rydberg Cs(nl) atom, can 
be described by a three-body model shown in Fig. [TJ The A + and e~ are the core and the valence 
electron of the Rydberg Cs(nZ) atom respectively and B is the ground state Cs(6s) atom. To be 
consistent with the above mentioned atomic part, r is the distance between the valence electron 
and the Cs + core, while r' is the distance between the valence electron and the Cs(6^) atom. And 
the two Cs atoms are separated by a large distance R. The total Hamiltonian can be obtained 
by taking into account of all the two-body interactions, namely, the molecular Hamiltonian of 
long-range Rydberg CS2 molecule can be written as 

#mol = #atom + V AB {R) + V Be {r') + T AB , (4) 

where Vg e (Vab) is the interaction potential between B and e~ (A + and B) 

= <5b) 

and Tab is the kinetic energy of the two Cs atoms. Here a a is the polarizability of a ground state 
Cs(6.s) atom and Eqs. © is valid for long range interactions. Similar to the atomic effective core 
potential V & t om (r), the short range behavior of V B e{ r ') is complicated consequently it is directly 
obtained by ab initio calculations. While the short range behavior of Vab(R) is n °t needed in the 
present work since the prime interest of the present work is the long range Rydberg Cs2 molecules. 

Using Born-Oppenheimer approximation, the molecular Schrodinger equation can be sepa- 
rated into two parts and the molecular potential energy curve can be obtained from the electronic 
Schrodinger equation 

(^atom + V Be (r , ))\|/ e/ (r ! R) =Ed(R)yp a {T,R) (6a) 

V mol (R)=E e i(R)+V A B(R) (6b) 

In this work, the electronic Schrodinger equation Eq. (l6al) is solved by perturbation theory. For 
ultralong Rydberg molecules the interaction potential Vg f? (/) has small matrix elements between 
Rydberg states of i/ a tom due to short range characteristics of Vg f? (r / ) and widely delocalized Ry- 
dberg electron density, which can also be seen from the results section. Consequently it is valid 



to set //atom as zero order Hamiltonian and the interaction potential Vg e as perturbation. Since Vs e 
only leads to couplings between the originally non-degenerate atomic states, the corresponding 
potential curve to the first order perturbation theory reads 

V mol (R) = V AB (R)+E el (R) = V AB (R) + (^atomlM'-Ol^atom} (7) 

Eq. dTJ) is the potential energy with respect to the corresponding dissociation limit of the two 
Cs(6s)+Cs(n/) atoms. The final numerical integrations for the perturbation term are calculated by 

(^toml^rOl^atom) = fjj WO*** (p)F /w (0, cp)r 2 sinQdrdQdq 

= 2nff V Be (r')Dl 2 nl (r)Y 2 (rc-arccos ^-^ cos9 ^ ;0 J r' 2 sin 0' dr' dQ' 

for different R with r = VR 2 + r J2 — 2Rr' cos 0' to properly treat the singular point of the interaction 
Vseif) at / = (which behave as — l/r'). 



C. Numerical Diagonalization of CS2 Molecular Vibrational Hamiltonian 

The vibrational states of Rydberg Cs2 molecules can be obtained by diagonalization of the 
molecular radial/vibrational Hamiltonian. Having the knowledge of the molecular potential energy 
V mo i(R) the molecular rovibrational Hamiltonian reads 

h 2 1 d d J 2 
HZi = T A B + V mol (R) = -^Y 2 te R2 M+ Vm ° m + W 2 (8) 

with the reduced mass /u = \mc s and the total angular momentum T 2 . The total rovibrational 
wavef unction can be written as ^mo^R) = M. v j(R)Yj 1 with M. v j(R) being the radial/vibrational 
wavefunction and Y" 1 the spherical harmonics as mentioned in previous section. The correspond- 
ing rovibrational energy level is E v j with vibrational quantum number v and rotational quantum 
number J. For a given angular momentum /, Eq. ([8]) can be diagonalized to give eigenenergy 
E v j and radial/vibrational wavefunction ~§t v j(R). Numerical diagonalization of Eq. © is per- 



formed in a Fourier grid discrete variable representation[|25ll3l|] to give vibrational energy E v q and 
corresponding vibrational wavefunction. 
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III. RESULTS AND DISCUSSION 



A. The Atomic Wavefunctions of Cs(ns/np) Rydberg States 

A typical Rydberg electron density is shown in Fig. [2] (for Cs(30p) state as an example). In the 
inset of Fig. |2] the Cs effective core potential VatomM is also shown. The effective core potential 
Vatom(f) is calculated by Gaussian09 suit of programs for r < 50 ao with a step size of 0.1 ao at 
CCSD(T)/QZVP level of theory. While for r > 50 a the V atom (r) is first calculated by Gaussian09 
at the same level of theory then fitted to Eq. ©. The fitted core polarizability for Cs + is a c = 
— 16.0 a,Q with a relative fitting error of less than 0.01%. The calculated core polarizability agree 
well with the experimental value of -15.5443 reported in Ref. J32]. For further calculations 
the experimental core polarizability a c = —15.5443 is aopted for r > 50 ao- The ab initio 
calculated effective core potential for r < 50 ao and the analytical one for r > 50 ao are connected 
smoothly as shown in the inset of Fig.|2l which generates the full effective core potential VatomM 
used for the calculations of atomic wavefunctions. 

The ns and np states of Cs atom wavefunctions are calculated by numerical integration of 
Eq. © with a step size of 0. 1 ao an d initial condition of the wavefunction at r > In 2 ao which 
obey asymptotic form 9t„/(r) °= }exp (—y/—2mE n ir/fi). As can be seen from Fig. [2] the Rydberg 
electron density is highly oscillating. The calculated most probable radius for the valence electron 
of the Cs(30p) state is 1297.7 ao which is the typical order of magnitude of the size of a Rydberg 
atom. For the Cs(30s), Cs(60s) and Cs(60p) states, the calculated most probable radii are 1250.9 
ao, 5989.4 ao and 6093.0 ao, respectively. The size of a Rydberg atom increases rapidly with 
increasing principle quantum number n and slightly increases with increasing angular momentum 
quantum number I, both due to increasing eigenenergy E n i. 



B. The Potential Energy Curves of Cs{ns / np)+Cs{6s) Rydberg Molecules 

Having the wavefunctions of Cs atomic Rydberg states at hand, the corresponding Rydberg 
CS2 molecular potentials V mo \(R) can be obtained according to Eq. © with the knowledge of the 
interaction potentials Vg e (/) and Vab(R)- The long range behaviors of Ve e (r / ) and Vab(R) can be 
well described by the ground state polarizability a a of Cs atom according to Eq. ©. The Cs atom 
polarizability a a has been obtained in the same way as the above mentioned core polarizability a c , 
namely, by fitting of the ab initio interaction energy Vs e (/) of r 1 > 50 ao calculated by Gaussian09 
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suit of programs at CCSD(T)/QZVP level of theory. The fitted ground state cesium atom polariz- 
ability is a a = —416.3 with a relative fitting error of less than 0.002%. The atom polarizability 
a a is significantly larger than the core polarizability a c since the outmost 6s valence electron of 
the Cs atom can be easily polarized. For r' < 50 ao, Vse(r') has been calculated by Gaussian09 
suit of programs at CCSD(T)/QZVP level of theory with a step size of 0.1 ao- 

The calculated potential energy curves of Rydberg Cs2 molecules are shown in Fig. [3] for 
Cs(30s /30p)+Cs(6s) and Fig.@]for Cs(60s/60p)+Cs(6s). The potential curves also oscillate due 
to the oscillating characteristics of the atomic Rydberg electron density. For details, the potential 
curves dissociating into 6s ground state and ns (np) Rydberg state are shown in the upper (bottom) 
panels of Fig. [3] and Fig.|4j Due to increasing size of a Rydberg atom with angular momentum 
quantum number /, each bottom panel (np+6s) has a larger equilibrium internuclear distance for 
the outmost potential well than that of the upper panel (ns+6s). In general the Rydberg molecules 
can be associated to the outmost potential well more efficiently than other potential wells. As far 
as the outmost potential well is concerned, the dissociation energy D e of the bottom panel (np+6s) 
is found to be considerably larger than the upper panel (ns+6s) which reflects stronger interaction 
between np and 6s states than the one between ns and 6s states for a given principle quantum 
number n. The calculated results for Cs*+Cs agree with the reported trends for Rb*+Rb in Ref. 
iHl. 

Details of the dissociation energies D e and equilibrium internuclear distance R m i„ for the out- 
most potential wells of Rydberg Cs2 molecules associated from different atomic states are shown 
in Table H Following the same symbols adopted throughout the present work, a Rydberg Cs2 
molecule is formed starting from a ground state Cs (6s) atom and a Rydberg state Cs (nl) atom. 
Detailed investigations are performed for the principle quantum number n varying from 30 to 70. 
As can be seen, the dissociation energy D e decreases almost in an exponential way with increasing 
principle quantum number n of the Rydberg atom which agrees with the experimentally observed 
trends for Rb*+Rb in Ref. Q4Q . This is because the Rydberg electron density becomes more de- 
localized as the principle quantum number n becomes larger. As a result the outmost maximum 
of the electron density, namely the most probable radial probability, decreases with increasing n 
which weakens the interaction between the two Cs atoms. While the molecular size, which is 
reflected by R m i n for the outmost potential well, increases according to the well known n 2 rule. By 
fitting the data reported in Table U one finds the relations R m j„ = l.65n 2 for ns+6s Rydberg Cs2 
molecules and R m i n = 1.68n 2 for np+6s ones, which are quite close to each other. 
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C. Vibrational States of Cs(ns/np)+Cs(6s) Rydberg molecules 

To further understand properties of the ultralong Rydberg Cs2 molecules, the molecular vibra- 
tional energy levels and wave functions are calculated based on the above mentioned potential 
energy curves. As mentioned in the method section, numerical diagonalization of Eq. © is per- 
formed in a Fourier grid discrete variable representation ll25ll to give vibrational energy E v q and 
corresponding vibrational wavefunction. Table |n] shows some vibrational energies E v q of the Ry- 
dberg CS2 molecule dissociating into Cs(30s) and Cs(6s) states. A total of 15925 grids equally 
spaced from R=107.5 ciq to R=1699.9 ao with a step of 0.1 ao, which covers 17 potential wells, 
are adopted for numerical diagonalization. The diagonalization finds 58 vibrational bound states 
among which 10 states are predominately located in the outmost potential well. Also shown in 
Table [II] are the average internuclear distance (R) and variance (AR) = a/ ({R— (R)) 2 ) of each 
corresponding vibrational state. In general the Rydberg molecules can be associated to the out- 
most potential well more efficiently than other potential wells. Therefore Table [TT] shows all the 
vibrational states which are predominately located in the outmost potential well. In this work all 
the vibrational energies are reported with respect to the corresponding dissociation energy. As can 
be seen the most deeply bound state in the outmost potential well is located at -498.5 MHz with 
an average internuclear distance (R) of 1246.6 ciq. Compared to the corresponding D e = —538.2 
MHz a relatively large zero point energy of about 40 MHz is found and the equilibrium vibrational 
frequency in this potential well is about 78 MHz according to the corresponding fundamental 
excitation energy. 

Similar calculations have been done for the Rydberg Cs2 molecules dissociating into Cs(30p) 
and Cs(6s) states. A total of 15251 grids equally spaced from R=174.8 ao to R=1699.8 ao with 
a step of 0.1 ao, which covers 15 potential wells, are adopted for numerical diagonalization. In 
total 80 vibrational bound states are found among which 17 ones are predominately located in the 
outmost potential well. Table [TIT1 shows all the 17 vibrational states in the outmost potential well. 
The most deeply bound vibrational states are found to be located at -1391.0 MHz with average 
internuclear distance (R) of 1292.5 ao- Compared to the corresponding D e = —1455.9 MHz the 
zero point energy is about 65 MHz and the equilibrium vibrational frequency in this potential well 
is about 128 MHz according to the corresponding fundamental excitation energy. 

In general, as the vibrational energy level increases the mean distance (R) and variance (AR) 
increase due to approaching of the dissociation limit as can be seen in Tables HH and Hill How- 
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ever, for the highest few vibrational states in Tables [TT] and [III they already have considerable 
population in the inner potential wells leading to quite large variances (AR) and decrease of the 
mean distances (R). Comparison between the above mentioned vibrational energies and the ones 
obtained from diagonalization of a relatively smaller Hamiltonian matrix which only covers the 
outmost potential well reveals that the vibrational levels in the outmost potential well can also 
be accurately obtained by the latter calculation (diagonalizing a smaller matrix). The root means 
square deviations between this two calculations are found to be less than 0.02 MHz for both the 
lowest 7 states in Table HI and the lowest 12 states in Table Hill If only the well bounded (lowest 
few) states are considered the deviations between the above mentioned two calculations are almost 
zero. To check the convergence of the above calculations, the vibrational energies in Tables HI and 
HIH are further compared with results from diagonalization of a larger Hamiltonian matrix which 
covers the long range limit up to R = 4000 ciq. The root mean square deviations between the two 
calculations are about 0.03 MHz for the vibrational energies in Table HI and 0.1 MHz for the ones 
in Table Hill Again the deviations are almost zero if only the well bounded (lowest few) states are 
considered. 

Based on the above mentioned error analysis, further investigations on vibrational states of Ry- 
dberg Cs2 molecules dissociating into Cs(6.s) and Cs(n/) with high principle quantum number n are 
performed by covering only the outmost potential well. For the Cs(60s /60p)+Cs(6s) molecules 
13583/12587 grids equally spaced from R=564 1.7/574 1.3 a to R=6999.9 a with a step of 0.1 a 
are adopted. Table HY] shows 3 vibrational bound states in the outmost well for Cs(605)+Cs(65) 
case and 5 vibrational bound states in the outmost well for Cs(60jc)+Cs(65) case. The corre- 
sponding zero point energy and equilibrium vibrational frequency in the outmost potential well of 
Cs(60s/60p)+Cs(6s) are 1.5/2.6 MHz and 2.7/4.9 MHz, respectively. 

IV. CONCLUSIONS 

Ultralong range Rydberg molecules have been extensively studied in the past decade due to its 
unique properties such as huge size, large dipole moment and long lifetime. In the present work 
ultralong Rydberg Cs2 molecules formed by Cs(ns / np)+Cs(6s) for (n = 30 — 70) have been inves- 
tigated focusing on molecular potential curves and vibrational properties. The Rydberg molecular 
potential curves are calculated by perturbation theory based on accurate Rydberg atom wave- 
functions. The Rydberg atom wavefunctions are generated from realistic ab initio effective core 
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potential rather than some empirical effective core potential, which serves as a good starting point 
for further calculations of potential energy curves. The calculated potential energy curves have 
oscillating nature due to the oscillation of the Rydberg electron density. The investigated Ryd- 
berg Cs2 molecular size increases with n according to R m i„ = 1.65n 2 for Cs(ns)+Cs(6s) case and 
Rmin = 1.68n 2 for Cs(np)+Cs(6s) case. The dissociation energy of the outmost potential well 
decreases almost exponentially with increasing principle quantum number n but increases dramat- 
ically with increasing angular momentum quantum number /. The largest binding energy of the 
vibrational ground states studied in this work is 1391 .0 MHz for Cs(30j!7)+Cs(65) case and quickly 
decreases to a few MHz for n = 10 case. According to our best knowledge this is the first time 
that vibrational properties of Cs*+Cs have been investigated quantitatively at the ab initio level. 
Consequently the reported results are expected to play important roles for future experiments on 
observation and manipulation of Cs*+Cs kind of molecules. For each vibrational levels predomi- 
nately located in the outmost potential well, apart from the vibrational energies, the mean values 
and variances of the internuclear distances are also reported which may provide further clues for 
the experimentalists. 
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Figure Legends 

Fig.® 

Rydberg molecule structure and coordinate system. The A + and e~ are the core and the valence 
electron of the Rydberg Cs(n/) atom respectively and B is the ground state Cs(6s) atom. See the 
text for details. 

Fig. m 

The radial probability density p(r) = |rlH„/| 2 of Cs(30p) atom. The effective core potential 
VatomM °f me Cs atom is shown in the inset. 

^•El- 
Molecular potential curves for Cs(305)+Cs(65) (Upper) and Cs(30p)+Cs(65) (Bottom). 

Fig.m 

Molecular potential curves for Cs(605)+Cs(65) (Upper) and Cs(60p)+Cs(6s) (Bottom). 
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TABLE I: The calculated dissociation energy D e and equilibrium internuclear distance /?„,„, for the outmost 
potential well of a Rydberg Cs2 molecules associated from a Rydberg Cs(ns/np) atom and a ground state 
Cs(6s) atom. 





ns + 6s 


np + 6s 


// 


R m in (P-q) 


D e (MHz) 


Rmin i. a 0) 


D e (MHz) 


30 


1244.6 


-538.2 


1291.4 


-1455.9 


35 


1787.0 


-197.3 


1843.2 


-542.7 


40 


2427.9 


-84.1 


2493.5 


-234.0 


45 


3167.4 


-40.0 


3242.5 


-112.5 


50 


4005.7 


-20.8 


4090.2 


-58.8 


55 


4942.7 


-11.5 


5036.8 


-32.8 


60 


5978.6 


-6.8 


6082.2 


-19.3 


65 


7113.5 


-4.2 


7226.5 


-11.9 


70 


8347.2 


-2.7 


8469.8 


-7.6 
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TABLE II: Vibrational properties of Cs(30s)+Cs(6.?) Rydberg molecules. E v q, (R) and (AR) are the vibra- 
tional energy, mean value and variance of the internuclear distance R, respectively. 



£ v0 (MHz) 


-498.5 


-420.5 


-347.4 


-279.6 


-217.3 


-160.9 


-111.2 


-69.2 


-32.7 


-10.1 


(R) (flo) 


1246.6 


1250.8 


1255.8 


1261.5 


1268.4 


1276.6 


1286.5 


1285.9 


1258.1 


1328.0 


(AR) (ao) 


18.6 


32.9 


43.8 


53.6 


63.3 


73.4 


85.3 


122.0 


215.5 


218.4 
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TABLE III: Vibrational properties of Cs(30/?)+Cs(6s) Rydberg molecules. E v $, (R) and (AR) are the vibra- 
tional energy, mean value and variance of the internuclear distance R, respectively. 



E v0 (MHz) 


-1391.0 


-1262.7 


-1138.9 


-1019.8 


-905.5 


-796.2 


-692.1 


-593.4 


-500.4 


(R) M 


1292.5 


1295.0 


1297.7 


1300.6 


1303.8 


1307.3 


1311.1 


1315.5 


1320.4 


(AR) (ao) 


14.5 


25.5 


33.4 


40.2 


46.5 


52.4 


58.3 


64.1 


70.1 


E vQ (MHz) 


-413.2 


-332.4 


-258.3 


-191.8 


-132.8 


-86.1 


-38.5 


-8.9 




(R) (ao) 


1325.9 


1332.3 


1339.2 


1345.3 


1351.0 


1256.0 


1268.8 


1247.9 




(AR) (a ) 


76.2 


82.8 


91.7 


101.3 


120.4 


225.3 


299.4 


377.1 
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TABLE IV: Vibrational properties of Cs(60s)+Cs(6s) and Cs(60/?)+Cs(6s) Rydberg molecules. E v0 , (R) 
and (AR) are the vibrational energy, mean value and variance of the internuclear distance R, respectively. 



60s + 6s 


60p + 6s 


E v0 (MHz) 


-5.3 


-2.6 


-0.7 


E v0 (MHz) 


-16.7 


-11.8 


-7.5 


-4.0 


-1.3 


(R) (a ) 


5995.5 


5930.9 


6076.7 


(R) (ao) 


6092.2 


6116.9 


6149.9 


6167.6 


6076.6 


(AR) (a ) 


100.3 


352.6 


457.4 


(AR) (ao) 


74.0 


133.7 


185.6 


272.2 


552.8 
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FIG. 1: Liu et al 
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FIG. 3: Liu et al 
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